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ABSTRACT 

Spatial intermittency in decaying kinetic Alfven wave turbulence is investi- 
gated to determine if it produces non Gaussian density fluctuations in the in- 
terstellar medium. Non Gaussian density fluctuations have been inferred from 
pulsar scintillation scaling. Kinetic Alfven wave turbulence characterizes density 
evolution in magnetic turbulence at scales near the ion gyroradius. It is shown 
that intense localized current filaments in the tail of an initial Gaussian probabil- 
ity distribution function possess a sheared magnetic field that strongly refracts 
the random kinetic Alfven waves responsible for turbulent decorrelation. The re- 
fraction localizes turbulence to the filament periphery, hence it avoids mixing by 
the turbulence. As the turbulence decays these long-lived filaments create a non 
Gaussian tail. A condition related to the shear of the filament field determines 
which fluctuations become coherent and which decay as random fluctuations. The 
refraction also creates coherent structures in electron density. These structures 
are not localized. Their spatial envelope maps into a probability distribution 
that decays as density to the power —3. The spatial envelope of density yields a 
Levy distribution in the density gradient. 

Subject headings: ISM: electron density — ISM: general — MHD — Turbulence 



1. Introduction 



Pulasr radio signals probe fluctuations in the local i nterstellar mediurn (|Armstrong et al. 



198ll ). The broad electron density fluctuation spectrum (lArmstrong et al.lll995l ) is commonly 
interpreted as a turbu lent inertial rang e . The pulsar signal width yields information about 
fluctuation statistics (IBhat et al.l |2004J : ISuttonI Il97ll ). The width scales as {R is the 
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distan ce to source), a result that is incompatible with Gaussian statistics flBoldyrev &: Gwinn 
2003af). The latter would p roduce a scaling of i?^, while R'^ is recovered for Levy statistics 
(iBoldyrev fc Gwinij l2003bl ) . A Levy-distributed random walk typically consists of a series 
of small random steps, punctuated by occasional Levy flights in which there is a single large 
jump to a new locale. In the context of a pulsar radio signal propagating through a Levy 
distribution of electron density fluctuations, a sea of low intensity density fluctuations would 
scatter the signal through a series of small angles. Intermittently, as the signal traversed an 
intense, localized density fluctuation, it would scatter through a much larger angle. 

The assertion that pulsar signals are dispersed by Levy-distributed fluctuations is a 
statistical ansatz validated to some degree by observation. This ansatz does not address 
the difficult and important question of what processes or conditions produce the statistics. 
It has been suggested that Levy st atistics can emerge from radio signal trajectories graz- 
ing the surface of molecular clouds (IBoldyrev fc Konigll 120061 ). Here we examine a different 
mechanism rooted in the turbulent cascade implied by the broad fluctuation spectrum. The 
mechanism is intrinsic spatial intermittency, a process known to create non Gaussian tails 
in the probability distribution function (PDF). In Navier-Stokes turbulence, intrinsic in- 
termittency takes the form o f randomly dispersed, localized vortex strands, surrounded by 
regions of relative inactivity ( Kerrlll985l ). Inte rmittency is most pr onounced at small scales. 
Intermittency also occurs in MHD turbulence (iGrappin et al.lll99ll ). However the statistical 
properties of electron density fluctuations in magnetic turbulence are not known. In this 
paper we address the fundamental and nontrivial question of whether electron density can 
become intermittent in the magnetic turbulence of the interstellar medium. The effect on 
pulse-width scaling requires that additional issues be addressed, and will be taken up later. 

The question of intermittency in pulsar scintillation is twofold. First, can intermit- 
tent electron density fluctuations in interstellar turbulence achieve the requisite intensity to 
change the PDF? To some extent this question has been answered by studies that show that 
passive advection and the limitations it places on electron density excitation (as indicated, 
for example, by mixing length arguments) applies only to scales larger than tens of gyro- 
radii. At smaller scales the electron density becomes active through kinetic Alfven wave 
(KAW) interactions with magnetic fluctuatio ns, exciting the internal energy to equipartition 
with the magnetic energy (ITerry et al.ll200ll ). Evidence for a transition to KA W dynamics 



near the gyroradius scale has recently been inferred from solar wind observations (IBale et al. 



20051 ) . Since scintillation is dominated by small scales, the regime of kinetic Alfven interac- 
tions is appropriate for studying the intermittency potentially associated with the scaling of 
the pulsar signal width. The second aspect of intermittency in the context of pulsar signals 
deals with how isolated structures can form against the homogenizing influence of turbulent 
mixing in a type of turbulence that does not involve flow. Virtually all mechanisms proposed 



- 3 - 



for intermittency involve flow or momentum, yet, the flow of ions in magnetic turbulence 
decouples from small-scale kinetic Alfven waves, with the interaction of magnetic field and 
density taking place against a background of unresponsive ions. 



While intermittency has been wid ely studied in hydrodynamic turbulence ( lFrischlll995l ) 
and MHD (jPolitano fc Pouquetlll995l ). historically the emphasis has been on structure and 
statistics, not mechanisms. Structure studies have included efforts to visualize intermittent 



structures (IHead fc Bandyopadhyaylll98ll ). Quantitatively, measurements of structure func 
tion scalings have been made to gauge how intermittency changes with scale (IShe fc Leveque 
1994 : iPolitano fc Pouquetlll995l ) . Statistical characterizations of intermittency generally pos- 
tulate a non Gaussian statistical ansatz, and the resultant properties are compared with 
measurements to determine if the ansatz is reasonable. These approaches do not address 
the mechanisms that endow certain fluctuation structures either with individual longevity or 
collective promi nence, in a sta tistical sense, relative to other regions in which such structures 
are not present (IWaleffelll997l ). The mechanistic approach is nascent but has already yielded 
signif icant insights in to the long-standing problem of subcritical instability in plane Pouseille 
flow jHof et ahlbnol . 



A starting point for our considerations are simulations of decaying KAW turbulence that 
showed t he emergence of coher ent, longlived current filaments under collisional dissipation of 
density (ICraddock et al.lll99ll ). In these simulations finite amphtude fluctuations in density 
and magnetic field decayed from initial Gaussian distributions. (The current, as curl of 
the magnetic field, was also Gaussian initially.) The distribution of current became highly 
non Gaussian as certain current fluctuations persisted in the decay long past the nominal 
turbulent correlation time. The longevity of these filaments enhanced the tail of the PDF, 
steadily increasing the value of the fourth order moment (kurtosis) significantly above its 
Gaussian value. While the PDF was affected by mutual interactions of filaments later in the 
simulation, initially the tail enhancement was dominated by the interaction of filaments with 
surrounding turbulence, and the lack of mixing of those filaments relative to the rapidly- 
decaying surrounding turbulence. Intermittency was not reported when resistivity dominated 
the dissipation. While these simulations showed intermittency in KAW turbulence, non 
Gaussian statistics was demonstrated for current fluctuations, not density. The turbulence 
decayed via collisional dissipation of density; the current had no direct damping. It is 
not clear what effect this had on density structure formation within the constraints of the 
resolution of the simulations. The question of intermittency in density therefore remains 
open. No mechanism for the intermittency was proposed. 



In this paper we will examine analytically the dynamics of structures in density and 
current and determine how one relates to the other. We will use analysis tools and results 
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developed t o understand the e mergence of long lived vortices in decaying 2D Navier-Stokes 
turbulence (iMcWilliamsl Il984j ) . For that problem, two time-scale analysis showed that the 
vortices are coherent and long lived becau se strong shear flow in the out er part of the vortex 
suppresses ambient mixing by turbulence (jTerrylll989l : iTerry et al.lll992l ). The ambient mix- 
ing would otherwise destroy the vortex in a turnover time. This mechanism for maintaining 
the coherent vortex in decaying tur bulence corr ectly predicts the observed distribution of 
Gaussian curvature of the flow field (jTerryll2000l ). 



We use two-time-scale analysis to describe coherent structure formation in decaying 
KAW turbulence. The following are obtained. 1) We identify the mechanism that allows 
certain current filaments to escape the turbulent mixing that otherwise typifies the turbu- 
lence. Current and density are mixed by the random interaction of kinetic Alfven waves. 
This process is disrupted in current filaments whose azimuthal field has an unusually large 
amount of transverse shear. This creates a strong refraction of turbulent kinetic Alfven waves 
that localizes them to the periphery of the filament and restricts their ability to mix current 
and density. 2) We derive a shear threshold criterion based on this mechanism. It identifies 
which current filaments escape mixing and become coherent, or long lived. The criterion re- 
lates to the Gaussian curvature of the magnetic field, providing a topological construct that 
maps the intermittency in a way analogous to the flow Gaussian curvature of decaying 2D 
Navier-Stokes turbulence. 3) We trace the relative effects of the refractive shear mechanism 
on current, magnetic field, density, and flux. The magnetic field and density have long-lived, 
localized fluctuation structures that coexist spatially with localized current filaments. How- 
ever, the magnetic field extends beyond the localized current. Like the magnetic field of a 
line current, it falls off as r^^. Because the density is equipartitioned with the magnetic 
field in KAW turbulence a similar mantle is expected for the density. This mantle tends to 
prevent the density kurtosis from rising to values greatly above 3; however, it is responsible 
for giving the PDF of density gradient a Levy distribution. 

This paper is organized as follows. Section II presents the kinetic Alfven wave model 
used in this paper. The two-time-scale analysis is introduced in Sec. III. Section IV derives 
the condition for strong refraction, and the resultant refractive boundary- layer structure 
for turbulent KAW activity in and around the coherent filament. The turbulent mixing 
stresses are determined in Sec. V, from which the filament and density lifetimes can be 
derived. Section VI discusses the Gaussian curvature and spatial properties of the current 
and density structures. The latter are used to infer heuristic PDFs. Conclusions are given 
in Sec. VII. 
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2. Kinetic Alfven Wave Model 

The shear-Alfven and kinetic-Alfven physics described in the introduction is intrinsic 
to models of MHD augmented by electron continuity. When there is a st rong mean field , 



the nonlinear MHD dynamics can be represented by a reduced description (jHazeltindll993l ). 
given by 



d'lb - CI 
^ + V||0 = r^J + Viin + t/— V^A X z • V^o, 

d 



^Vi0 - V0 X z ■ VVi0 = -V|| J, (2) 

dfi ^ - C 1 - 

— - V0 X z- Vn + ViiJ- V0 X z- Vno = 0, (3) 

at " VAriQ 

where 

d 

V|| = — + VV^ X z ■ V, (4) 

and J = = d'^ip/dx'^ + d'^ip/dy'^. In the reduced description the perturbed magnetic 

field is perpendicular to the mean field and can be written as b/B = Vip x z, where z 
is the direction of the mean field, and = {Cg/ c)eAz/Tf. is the normalized parallel com- 
ponent of the vector potential. The fiow has zero mean and is also perpendicular to the 
mean field B. It can be expressed in terms of the electrostatic potential as — V0 x z where 
= (Cs/VA)e0/Te is the normalized electrostatic potential. The normalized density fluc- 
tuation is n = (Cs/V^)n/no, where no is the mean density, and rj = {c^ / 4:7!-VaPs)Vsp is the 
normalized resistivity, where rjgp is the Spitzer resistivity. Spatial scales are normalized to 
Ps = Cs/^i, time is normalized to the Alfven time ta = Ps/Va, Cg = {Te/rriiY^'^ is the ion 
acoustic velocity, Va = B / {AnminQy/'^ is the Alfven velocity, and VLi = eB/rriiC is the ion 
gyrofrequency. Within their limitations (isothermal, incompressible fluctuations), Eqs. (1)- 
(4) are valid for scales both large and small compared to the gyroradius, as well as the 
intermediate region. 

Equation (3) is the electron continuity equation. The advective nonlinearity, V0xz- Vn, 
couples electron density fluctuations to the flow. If there is a nonuniform mean density, ad- 
vection drives weak density fluctuations of amplitude h ^ (5/L„)no, where 6 is the scale 
of density fluctuations and L„ is the mean density scale length. The continuity equation 
also contains a compressible nonlinearity, V|| J, whereby compressible electron motion along 
magnetic fleld perturbations provides coupling to the magnetic fleld. Electrons act on the 
magnetic fleld through parallel electron pressure in Ohm's law, expressed as V||n in Eq. (1). 
The couplings of magnetic fleld and density are weak at scales appreciably larger than the 
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ion gyroradius. On those scales the advection of electron density is passive to a good ap- 
proximation, and governs electron density evolution. In th e region around 6 = lOps, the two 



nonlinearities in each of Eqs. (l)-(3) become comparable (iTerry et al.ll200ll ). For 6 < lOps, 
V||n begins to dominate V||0 in Eq. (1), and V|| J begins to dominate Vcf) x z- Vn in Eq. (3). 
This is a very different regime from incompressible MHD, where the magnetic field and flow 
actively exchange energy through shear Alfven waves. In a turbulent cascade approaching the 
ion gyroradius scale from larger scales, the energy exchanged between flow and magnetic field 
in shear Alfven interactions diminishes relative to the energy exchanged between the electron 
density and the magnetic field through the compressible coupling. Consequently flow decou- 
ples from the magnetic field, increasingly evolving as a go-it-alone Kolmogorov cascade, while 
electron density and magnetic field, interacting compressively through kinetic Alfven waves, 
supplant the shear Alfven waves. Once the kinetic Alfven wave coupling reaches prominence, 
the internal and magnetic energies become equipartitioned, J n^dV ~ j \W'4)\'^dV , even if 
the internal energy is only a fraction of the magnetic energy at larger scales. If there is no 
significant damping at the ion gyroradius scale, the large-scale shear Alfven cascade contin- 
ues to gyroradius scales and beyond though kinetic Alfven waves. The gyroradius scale at 
which KAW dynamics is active is order 10* cm in the warm ISM. This is small relative to the 
scale of i ntermittent flow struct ures in molecular gas clouds, recently reported to be order 
10^* cm (IHily-Blant et al.l 120071 ). This scale difference is crudely consistent with the high 



magnetic Prandtl number of the warm ISM. The value Pr ~ 10^^ allows very small scales 
in the ionized medium, before dissipation becomes important, relative to scales of viscous 
dissipation in the clouds. The gyroradius scale of intermittent KAW structures makes direct 
visualization in the ISM difficult. 

In the KAW regime, the model can be further simplified by dropping the flow evolution. 
This leaves a KAW model in which electron density and magnetic field interact against a 
neutralizing background of unresponsive ions. 



dtb dn ^ ^ ^ c Cs I „ 7 „ 

^ = + V^p X z-Vn + r]J + — Vip x z ■ Vkq, 

ot oz Va riQ 



— + — + VV^xz-VJ- uV^n 
at oz 



0. 



(5) 
(6) 



Solutions of this model cl osely approximate those of Eqs. (l)-(3) when the scales are near 
the gyroradius or smaller (ITerry et al.lll998l ). This model assumes isothermal fluctuations, 
consistent with strong parallel thermal conductiv ity. Equations (5) and (6) are fluid equa- 
tions, hence Landau- resonant (IHowes et al.ll2006l ) and gyro-resonant dissipation, which may 
be important in the ISM, are not modeled. Ohm's law has resistive dissipation, and density 
evolution has collisional diffusion. Depending on the ratio rj/fi, either of these dissipation 
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mechanisms can damp the energy in decaying turbulence, however, the damping occurs at 
small dissipative scales. We will focus on inertial behavior at larger scales. We assume that 
mean density is nearly uniform, and neglect the last term of Eq. (5). 

The dispersion relation for ideal kinetic Alfven waves is determined by linearizing 
Eqs. (5) and (6), neglecting resistive dissipation r^J, and introducing a Fourier transform 
in space and time. The result is = kzk±, where k_|_ ■ z = 0, or ± k^z. If dimen- 
sional frequency and wavenumbers oj, k^, and k± are reintroduced, the dispersion relation 
is a) = VAkzk±ps. The wave is seen to combine Alfvenic propagation with perpendicular 
motion associated with the gyroradius scale. The KAW eigenvector yields equal amplitudes 
of magnetic field and the density, —k±ipk = ibk = rik, with a phase difference of 7i/2. 

In magnetic turbulence with its hierarchy of scales, kinetic Alfven waves also propagate 
along components of the turbulent magnetic field. In the reduced description the turbu- 
lent field is perpendicular to the mean field, hence the dispersion relation of these kinetic 
Alfven waves carries no kz dependence. To illustrate, we isolate such a fluctuation from the 
mean-field kinetic Alfven wave by setting k^ = 0; with this wavenumber zero, we drop the 
subscript from k^; we consider a turbulent magnetic field component h^^jBo = ikg x zip^^^ 
at wavenumber ko that dominates the low-fc fluctuation spectrum; and we look at the dis- 
persion for smaller scale fluctuations satisfying k ^ ko. The latter conditions linearize the 
problem, yielding a dispersion relation for kinetic Alfven waves propagating along the tur- 
bulent field b/jQ according to u; = iko x z ■ k ipkok = [b^^, ■ k/i?] k. Reintroducing dimensions, 
= Ki(bfco-k/i?)A;ps. We see that the dispersion is Alfvenic, but with respect to a perturbed 
field component that is perpendicular to the mean field. Hence the frequency goes like k'^ps 
instead of kzkps- 



3. Two Time Scale Analysis 

To understand and quantify the conditions under which a coherent current fluctuation 
persists for long times relative to typical fluctuations, we examine the interaction of the co- 
herent structure with surrounding turbulence and derive its lifetime under turbulent mixing. 
The interaction is describ ed using a two-tim e-scale analysis, allowing evolution on disparate 



time scales to be tracked (ITerry et al.lll992l ). The coherent structure, a current filament with 
accompanying magnetic field and electron density fluctuations, evolves on the slow time scale 
under the rapid-scale-averaged effect of turbulence. On the rapid scale the filament is essen- 
tially stationary, creating an inhomogeneous background for the rapidly evolving turbulence. 
Identifying conditions that support longevity justifies the two time scale approximation a 
posteriori. Simulations suggest the filament is roughly circular. If coordinates are chosen 
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with the origin at the center of the filament, a circular filament is azimuthally symmetric, 
while the turbulence breaks that symmetry. 

The filament current is locahzed, hence its current density becomes zero at some dis- 
tance from the origin. The localized current profile necessarily creates a magnetic field that 
is strongly inhomogeneous. On the rapid time scale over which the turbulence evolves, this 
field, which is part of the coherent structure, is essentially stationary. It acts as a secondary 
equilibrium field in addition to the primary equilibrium field (which is homogeneous and di- 
rected along the z-axis). Turbulence, in the form of random kinetic Alfven waves, propagates 
along both the primary and secondary fields. Because the primary field is homogeneous its 
effect on the turbulence is uninteresting. However, the secondary field is strongly sheared 
because of the local inhomogeneity created by the structure. Strong shear refracts the tur- 
bulent kinetic Alfven waves. In the subsequent analysis we will ignore the primary KAW 
propagation, which we can do by setting kz — 0, and focus on the refraction of KAW prop- 
agation by the secondary magnetic field shear. Strong refraction will be shown to localize 
kinetic Alfven waves away from the heart of the filament, allowing it to escape mixing and 
thereby acquire the longevity to make it coherent. 

With J = V^^V', we apply the separation of long and short time scales to fluctuations 
in n and as follows: 



where the symbol F represents either ijj or n, with ipo and Uq the flux function and density of 
the coherent structure, and i/j and n the turbulent fields of flux and density. The variables for 
slowly and rapidly evolving time are r and t. The origin of a polar coordinate system with 
radial and angle variables r and 6 is placed at the center of the structure. The structure 
is assumed to be azimuthally symmetric. The turbulence evolves in the presence of the 
structure, hence it is necessary to specify the radial profile of tpo, or more explicitly, the 
profile of the secondary, structure held V^/'o(r) xz — Be{r)9. As a generic proflle for localized 
current we adopt a reference proflle that peaks at the origin and falls monotonically to zero 
over a flnite radius a. For simplicity we take the variation as quadratic, giving 



Fo(r,T)+F(r,^,i), 



(7) 




(8) 



where Jo('^) — 



Integrating the current we obtain 




(r < a). 



(9) 



Be{r) 



4r 



(r > a). 



(10) 
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These profiles are for reference. Shortly we will introduce a more general description for 
a filament whose current peaks at the origin and decays monotonically. The current of 
the coherent filament is wholly localized within r = a. However, the magnetic field is not 
localized, but slowly decays as r^^ outside the filament. The quantities in Eqs. (8)-(10) 
all evolve on the slow time scale r. The dependence on r is not notatcd because when Bg 
appears in the turbulence equations, it is a quasi equilibrium quantity on the rapid time 
scale. 

To describe the rapid time scale evolution and its azimuthal variations we introduce a 
Fourier-Laplace transform. 



27ri 



+ioo+7o 
■joo+70 



<i7 X] Frnni^) exp[-im6'] exp[7t]. 



where 70 is the shift of the complex integration path of the inverse Laplace transform. The 
radial variation of BQ(r) creates an inhomogeneous background field for the turbulence, 
making Fourier transformation unsuitable for the radial variable. The Laplace transform 
is appropriate for turbulence that decays from an initial state. To obtain equations for 
the slowly evolving fields ■00 and uq, we average the full equations over the rapid time 
scale t. This is accomplished by applying the Laplace transform to the equations, and 
integrating over t. The integral selects 7 = as the time average, i.e.., J dtf{T,t) — 
J dt J d7(27ri)~^/y(r) exp(7t) = /y=o(T). Applying this procedure, the evolution equations 
for the slowly evolving fields are given by 



dr 



1 
27ri 



-ioo+70 



ioo+70 



d_ 

' dr 



V r 



nm',y ) = 0, (12) 



dn^=o{T) _^ J. 
dr 2m 



/—ioo+70 



' -ioo+70 



d 



br{-m',-Y)— 1" be(-m'.-'y') I 



dr 



X 



Id/ d \ 
r dr\ dr J 



m 



m ,7' 



0, 



(13) 



where 6j.(_^'__y), &e(_^'__y), n^',y, and ■0m' ,7' a-re understood to depend on the radial vari- 
able r; br(-m',-Y) — (~^"^V0V'-m',-7'; and bo(^-.rn' -y) — ~{d/dr)ip^m'-Y- The correlations 
{brdh/dr), (boh), {brdV'^ip/dr), and (beV'^ih), which appear in Eqs. (12) and (13), are turbu- 
lent stresses associated with random kinetic Alfven wave refraction. Their fast time averages 
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govern the mixing (transport) of the coherent fields. These stresses must be evaluated from 
solutions of the fast time scale equations to find the lifetime of the structure. 

The evolution equations for the rapidly evolving turbulent fiuctuations are 

fioo+7() 



r 



d i{m — m') dipra',') 



dr 



dr 



m— m ,7— 7' 



X 



— ^m,7^^o(r), (14) 



7n. 



■m.7 



^ \ r ) r dr\ dr ) 

18/ iJ \ 
r dr\ dr J 



2 /•joo+70 



X 



im' ~ d i{m — m') dtjjm' 
r dr r dr 



ioo+70 

r\2 



1 d / d \ {m — m') 



m—m ,7—7 



im~ d 

— V'm,7^^o(r). 



(15) 



We have not shown the dissipative terms in accordance with our focus on inertial scales. 
The sources contain gradients of no(r) and Jo{r). These arc the density and current of 
the coherent structure, but, unlike ?/''7=o(''") and n^=o(''"), are not evaluated in the Laplace 
transform domain. Three terms drive the evolution of dip/dt and dn/dt in each of these 
equations. The first term describes linear kinetic Alfven wave propagation along the inho- 
mogeneous secondary magnetic field Be of the coherent structure. The second term is the 
nonlinearity, and describes turbulence of random kinetic Alfven waves. The third term is 
proportional to mean-field gradients. It is a fluctuation source via the magnetic analog of 
advection (V0 X z • V = v • V — > V'lp x z ■ V). It yields quasilinear diffusivities for the 
turbulent mixing process. For example, if the kinetic Alfven wave and nonlinear terms of 
Eq. (15) are dropped, the solution is 



n. 



(i) 
m,7 



lim ~ dJoir) 
7 r dr 



(16) 



The superscript (i) indicates that, for deriving diffusivities, this density is to be substituted 
iteratively into the correlations of the turbulent stresses. From Eq. (12) these correlations 
are (prfi) = (br^n) + {brU^^^) and {bedn/dr) = (b^Q^dh/dr) + {bgdn'^^y dr) . Substitution of 
Eq. (16) yields mean turbulent diffusivities for ■0o- Similarly, if Eq. (14) is solved by dropping 
its kinetic Alfven wave and nonlinear terms, we obtain 



ii) 
m,7 



lim ~ dno{r) 
7 r dr 



(17) 
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Substituting this solution into the correlations {KV^ip) = {WV^'ijj)+{brV^tf(^) and (bedV'^ip/dr) = 
{hfdV^i)/dr) + {hedV'^^^^/dr) of Eq. (13), mean turbulent diffusivitics arc obtained for n^. 
Off-diagonal transport (relaxation of by gradient of no) can also be obtained by sub- 
stituting Eq. (17) into {hru) and {hgdn/dr). The role of the nonlinear and kinetic Alfven 
wave terms omitted from Eqs. (16) and (17) is to modify the time scale 7 and couple the 
sources. This is calculated in the next section. The inverse of 7 represents the lifetime of 
the correlations (6r^), (bedh/dr), (brV^ip), and (bedV'^'ip/dr). Generally the nonlinear terms 
enhance decorrelation, increasing the effective value of 7. If the shear in is strong, the 
kinetic Alfven wave term increases 7 even further. 

The role of shear in the kinetic Alfven wave terms is not explicit but should be, so that 
it can be varied independently of the field amplitude Bg{ro) at some radial location tq. In 
explicitly displaying the role of shear we note that if Bg{r) ~ r, as would be true if the 
current density Jq were uniform, the kinetic Afiven wave term is independent of r. In this 
situation the phase fronts of kinetic Alfven waves propagating along Bo are straight-line rays 
extending from the origin. Shear in Bg, occurring through nonuniformity of Jq, distorts the 
phase fronts, as shown in Fig. 1. Distortion occurs if Bg has a variation that is not linear. 
From Eq. (9) we note that the variation of Bg for our chosen structure profile is linear for 
r <^ a, with variations developing as r — >■ a. Therefore, it makes sense to quantify the shear 
by expanding Bg{r)/r in a Taylor series about some point of interest. Obviously, the shear 
is zero at the origin, and becomes sizable as r — > a. Expanding about a reference point Tq 
away from the origin. 



Bg{r) Bg{ro) d f Bg 

= + {r-ro)—[ — 

r tq ar \ r 



+ 



(18) 



If Bg{r) varies smoothly, as is the case for a monotonically decreasing current profile, we can 
truncate the expansion as indicated in Eq. (18) and use that expression as a general current 
profile. Looking at the kinetic Alfven terms of Eqs. (14) and (15), the first term will produce 
a uniform frequency that Doppler shifts 7 by the amount imBg{r-Q) /tq. The second term 
will describe KAW propagation in an inhomogeneous medium with its attendant refraction. 



4. Refraction Boundary Layer 

We rewrite Eq. (15), substituting the expansion of Eq. (18), yielding 



771, 



m,7 



im[r 



-)*(-)L/»'*-^2^/._/^'? 



X 



- 12 - 



im' ~ d 
r or 



i{m — m') dipr, 



dr 



V 



2 

m—m 



m—m' ,7—7 



Iffl ~ Q 

r or 



where 



1 d r d \ rn^ 
r dr\ dr / 



(19) 



(20) 



and 7 = 7 + imBg{rQ)/rQ. When d/dr[Bg/r]\rQ is large, the shear in Bg refracts turbulent 
KAW activity. The process can be described using asymptotic analysis. In the limit that 
d/dr[Bg/r]\rQ becomes large asymptotically, the higher derivative nonlinear term is unable 
to remain in the dominant asymptotic balance unless the solution develops a small scale 
boundary layer structure. The layer is a singular structure. Its width must become smaller 
as d/dr[B0/r]\ro becomes larger, otherwise the highest order derivative drops out of the 
balance and the equation changes order. This is the only viable asymptotic balance for 



d/dr[Be/r] 



00. The boundary layer width Ar is readily estimated from dimensional 



analysis by noting that r — ro ~ Ar, dnm{t)/dr ~ rim/^r, and treating d / dr[Bg / r]\ro = j' 
as the diverging asymptotic parameter. The asymptotic balance is 



a Ar 



if - 00), 



(21) 



yielding 



Ar 




(/ - 00). 



(22) 



The length Ar is the scale of fluctuation variation within the coherent current filament. In the 
simulations, the filaments were identified as regions of strong, localized, symmetric current 
surrounded by turbulent fluctuations. Consequently, Ar represents a fluctuation penetration 
depth into the structure. We derived the layer width Ar from linear and nonlinear kinetic 
Alfven wave terms operating on flux in the density equation. Identical operators apply to n in 
the flux equation. Hence Ar is the width of a single layer pertaining to both the density and 
current fluctuations of refracted KAW turbulence. This structure is shown schematically in 
Fig. 2. The above analysis indicates a single layer width and gives its value. It does not give 
the functional variation of current and density fluctuations within the layer, either relative 
or absolute. 

In the simpler case of intermittency in decaying 2D Navier-Stokes turbulence, statistical 
closure theory was used to derive spati al functions describ ing the inhomogeneity of turbulence 
in the presence of a coherent vortex (ITerry et al.lll992l ). There, coherent vortices suppress 
turbulent penetration via strong shear flow, analogous to role of refraction here. For the 
KAW system the closure equations are much more complicated and not amenable to the 
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WKB analysis that gave the functional form of the boundary layer in the Navier-Stokes 
case. However, the closure remains useful. It provides a mathematical platform from which 
to calculate all aspects of the interaction of filament and turbulence, including the accelerated 
decay of turbulence within the boundary layer, the spatial characteristics of the layer, and 
the amplitudes of n and ip. These are necessary for calculating turbulent mixing rates of the 
filament current and density. 

Closures can be applied to intermittent turbulence even though they rely on Gaussian 
statistics. The filaments, which make the system non Gaussian as a whole, arc quasi sta- 
tionary on the short time scale of turbulent evolution. Therefore, on that scale their only 
effect is to make the turbulence inhomogeneous. The fast time scale statistics are a property 
of fast time scale nonlinearity, and remain Gaussian. The closure equations are 

~ d'^ ~ d'^ ivn ~ d 

(23) 

^fira,^ - im{r - ro)j'V^^^,^ - D^^l (m, 7) ^ V^^^,^ - D^^l (m, 7) ^ V^^„,^ 

- 7)^^'V7 - ^)^^rnn = ^"^0 (r) . (24) 

This system is complex. The six diffusivitics all contribute to the lowest order as j' — > oo. 
(The diffusion coefficients and derivatives arc not of the same order, but their product is.) 
Moreover there is varied dependence on ffuctuation correlations, and there are complex 
turbulent decorrelation functions. For example 



where 



p — 

m. — 



-j^ /-too+io m' { / ~ m! 

+ Ki{m -m\-f- 7')P^i^,AW^,y (|V-m' -7'^m',7' ) |, (25) 



+ 1-D^lirn,^)§r2^l-D^lirnn)^], (26) 



-l + D.p4m,-f)—Vl-——no) [imj'{r-ro)+D^^{m,-f)—y (27) 

and Al^(7, 7') is the decorrelation rate for fluctuations at 7' driving 7. Expressions for the 
other diffusivities are given in the appendix. We now calculate turbulence properties from 
the closure equations. 
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4.1. Decay Time for Turbulence in the Filament 

The time scale of turbulent evolution in the filament is given by 7. As j' — > 00 this is 
dominated by the refraction. Hence the first terms of Eq. (23) and Eq. (24) must balance the 
second terms, which in turn, must balance the remaining nonlinear terms. If Eqs. (23) and 
(24) are solved jointly retaining the first two terms, 7 ~ im{r — ro)j'V as j' 00. Because 
(r — To) ~ Ar and V ~ 1/Ar, 

7 ~ imj' {j' — >• 00). (28) 

This time scale is purely imaginary, i.e., oscillatory, when derived from a balance with only 
the linear Alfven terms. When the diffusivities are included, it is complex. This rapid decay 
suppresses turbulence in the filament relative to levels outside the filament. 



4.2. Alfvenic Boundary Layer Width 

The width Ar, as derived previously, comes from independent balances in the equations 
for n and ip, and does not account for the kinetic Alfven wave dynamics that links n and ip. 
To do that, Eqs. (23) and (24) are combined into a single equation by operating on Eq. (24) 
with —im{r—rQ)i' — D^nijn, ^)d'^ jdr^ and substituting from Eq. (23). The resulting equation 
is eighth order in the radial derivative, and unsuitable for WKB analysis. However, we can 
determine the radial scale as j' 00 by dimensional anlaysis, taking ~ d"^ /dr"^ ~ 1/Ar^ 
and solving algebraically. This is the same procedure used to obtain Eq. fl22l) . Formally 
treating Ar as a small parameter, we account for the fact that the diffusion coefficients 
have different scalings with respect to Ar, based on different numbers of radial derivatives 
operating on quantities within the coefficients. Arbitrarily taking Dnn as a reference diffusion 
coefficient, the definitions in the appendix show that if we define D^^^ = Ard^^^, D^^ = 
{Ar)^S^^, D^nn = {ArySnl, d'H = Snl, D^^ = (Ar)2(i^^, and D^^ = Ard^n, then the lower 
case diffusivities d^^^, rf^^^, dnl, dnn, d^^, and d^n are all of the same order. We formally 
order the large parameter j' by taking mj' —>■ mj'/e and 7 — > 7/e, where the controlling 
asymptotic limit becomes e — > 0. The relationship between e and Ar will be derived by 
requiring that the asymptotic balance be consistent. After all leading order expressions are 
derived, e is set equal to 1. 

Substituting these relations into Eqs. (23) and (24) and solving, we obtain: 

^,,("^^ + ^(2)^ -(i, + _ j!_r^™,-'/'w« + c;(2)^ \ /, (1)^(2) 

"1/"/' y^nn ^ ^nn j ""pn ^ %^ ^ J [ {"'nip ^ ""nip ^ ""^n J 1 \^nn ^ "-nn 



Ar^ 
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—d 



+ m J +7 



m,7 



7 - 



d^^^ + d 



(2)n 

nn 



5,/ 



ArS„ 



where 



im 7 



d 



To 

im ~ d 
To dr 



no(r). 



(29) 

(30) 
(31) 



are the turbulence sources described in the previous section. The left hand side is a dimen- 
sional representation of a Green function operator that governs the response to the sources. 
The spatial response decays inward from the edge of the filament where both the sources 
and the shear in Bg are strong. Consequently, the field ipm,^ appearing in the sources Sn and 
is understood to be characteristic of the filament edge, and therefore ambient turbulence, 
while "^m,"/ is a response accounting for the the refractive decay inside the filament. The scale 
length of the response Ar is found by solving the homogeneous problem, i.e., by setting the 
left hand side equal to zero and solving for Ar. In the limit e ^ 0, turbulence remains in 
the dynamics and contributes to Ar only if e ~ Ar. Otherwise, the dynamics is laminar. 
The solution for Ar is 



Ar) ~ 2dl 2dl 



imj'd2 - 74) ' - 4(i? (m^f + 7') 



1/2 



(32) 



d^}l + (iiw, + d^n, and 4 



u-nn ~r u.nn <-l'ipip ■ 



where d-^ - d^^{dnn +dnn)- d^n{d^^ -r u,„^y , ^-2 - u,„^ -r u,„^ 

This is the Alfvenic generalization of Eq. (22). It is more complicated but gives identical 
scaling. Recalling that all the lower case diffusivities have the same scaling and replacing 
them with a generic d, the solution scales as e^/Ar^ ~ mj'/d. Setting e — 1, 



Ar 



mj' 



(33) 



The generic diffusivity d can be evaluated from the definitions given for specific diffusivities 
in the appendix. If the turbulent decorrelation functions are evaluated in a strong turbulence 
regime (turbulence time scales << linear time scales), d ~ ipm/O', reproducing Eq. (22). 



4.3. Boundary Layer Structure of Turbulence 



Although the structure function has not been solved (just its radial scale), its form in 
simpler cases illustrates the rapid decay of turbulence across the boundary layer, from the 
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edge inward. Where WKB analysis is possible, the leading order spatial Green function has 
the form 



G{r\r') - expl - C 



r< - rp 
\J d/mj' 



exp 



r> - rp 
V d/mj' 



(34) 



where C is a complex constant with positive real part, r< (r>) is the smaller (larger) of r 
and r', and a is a positive constant determined by the order of the homogeneous operator. 
Here our dimensional solution of the problem, carried out by inverting Eq. (29), captures 
the radial integral over a structure function like that of Eq. (34). Solving Eq. (29) we obtain 



ii) 
m,7 



7 



7 



+ imj + d. 



ArSr, 



m 



7 



V^(ro)K + ArJ^]. (35) 



The temporal and spatial response to turbulent sources and Sn at a point tq in the filament 
edge appears here as a structure factor of magnitude 7"^ multiplying the source. The product 
of source and response yields the value of tp^^^ inside the boundary layer. Beyond Ar the 
response decays with an envelope hke that of Eq. (34). The part of the source proportional 
to Uq is essentially larger than the part proportional to Jq by 0(a/Ar). However if ip^'^^ 
is substituted into the correlations of the equation for ■0o [Eq. (12)], the Jg-part yields the 
diagonal terms. The density is given by the dimensional representation of Eq. (23), 



imj'Ar + 



Ar 



-1 



D 



7V'' 



m,7 



V'm,7(?^o) n'o + Ar Jq 



aj'Ar 



(36) 



4.4. Condition for Strong Refraction 

The layer width Ar is both the embodiment of the strong refraction of turbulent KAW 
activity in the filament by the large magnetic field shear j', and a condition for the refraction 
to be sufficiently strong to modify the scales of turbulence in the filament relative to those 
outside it. With a the scale of typical fluctuations of interest, the refraction is strong when 
Ar/a <^ 1, or 

- ~ ,/4^ ~ \Ph « 1- (37) 
a y a^mf y a-^mj' 

As a condition for strong refraction it makes sense to use values for d or ip that are typical of 
the turbulence in regions r > a where there are no intense fllaments. Inside a strong fllament 
the reduction of turbulent KAW activity represented by the structure factor 7"^ makes 
d\r<a <^ d\r>a- Accordingly, the boundary layer width Ar/a is smaller than [d\r>a/ ci^iTij']^/^ . 
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5. Filament Decay from Mixing Stresses 



The long time evolution of the filament fields ipoir) and no(r) is governed by the mixing 
stresses of Eqs. (12) and (13). These can now be evaluated using the boundary layer responses 
ijjm,-y and nm,7 derived in the previous section. Because these fields are confined to the 
layer, the time scale r is a mixing time across the boundary layer. For the diagonal stress 
components, the mixing is diffusive. The asymptotic behavior of the boundary layer yields 
the following dimensional equivalents: r — ro ~ Ar and d/dr ~ 1/Ar, as before; J dr ^ Ar; 
J d'j r.-' mj'; i>'„i.,7 l^Jrn{t) /mj'\t=Q; and J d'y'{iljh) ~ '4'{t)n{t)/mj'\t=o. (The latter two 
expressions are inverse Laplace transform relations.) With these conventions 



9V'7=o(t) 1pj=o 



-1 /■-ioo+70 r • / , 



dr 



Tip 



+ 



dr 



d ^ 



'),'^m',7' 



(- 

\dr 



-7'^m',7' 



{be 



i=0 



\t=0 



rin 



ArJ^ 



(38) 



where bom — i^m/ ^f- The factor in the right-most form makes large, i.e., mixing 
across the boundary layer is impeded by refraction. The turbulent fields in these expressions 
are filament edge fields, i.e., they are characteristic of ambient turbulence. The mixing time 
for current can be obtained by operating with on both sides of the Eq. (38). On the left 
hand side V^'0-y=o 1^^=0/0^, while on the right hand side, 1/Ar^. Consequently, 



7=0 



E 



1 



{be 



i)\t=0 + (&em)lt=0 



Tin 



(39) 



This time scale is much shorter because current, as a second derivative of ^0, has finer scale 
structure. If the filament is Alfvcnic, i.e., tiq ^ aJo, the mixing time is dominated by the 
part of Eq. (39) that is proportional to Uq. This represents off-diagonal transport of current 
driven by density gradient. The diagonal transport (driven by Jg) is current diffusion, and is 
slower by a factor a/ Ar. In the discussions that follow we will deal with the current diffusion 
time scale, although similar behavior will hold for the off diagonal transport. The mixing 
time for density is 



9n^=o(r) n^=o 



dr 



E 



"Ia^^-^' 



(40) 



Here the dominant component (proportional to n'o) is diffusive. 



We evaluate these boundary layer mixing times relative to the two turbulent time scales 
of the system. These are 7"^, the turbulent decay time in the layer, and ta — a^/&6»|r->o) a 
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turbulent Alfven time outside the filament. To simplify expressions we note that Alfvenic 
equipartition implies that {hefi) ~ (6^). We also note that ^^=o is in the Laplace transform 
domain, whereas and Jq are in the time domain. Under the inverse Laplace transform, 
V''7=o7 ~ ''Pii, t) = i^o- The scale of the filament is a, so ip-y=ol ~ a'^Jo- Similarly, u'q ~ no/a, 
Jq ~ Jo/a, and j' = d/dr{Bg/r) Bg/a^. We assume the filament is Alfvenic, making 
Hq ~ gJq. With these relations, 

ArBl /a \3 /5fl\3/2 



= = _^.^_J .^^J . (41) 

The last two equahties make use of Eq. (37), and the fact that Ar/a and the mixing fluc- 
tuations are referenced to ambient turbulence levels for which bg ~ ip/a ~ br- Equation 
(41) indicates that turbulent diffusion times across the mixing layer Ar for both no and Jq 
are comparable, and arc much longer than the decay times of turbulence in the layer. The 
strong shear limit, previously indicated by j' oo, is here replaced by Bg — >■ oo, because 
with a fixed radius a, strong shear means large Bg. In terms of ta, 

ArBg / a \ (Bg\V^ 



a bg ^^rJ \b, 

indicating that these diffusion times are longer than the Alfvenic time of the ambient tur- 
bulence. 

Either of the above expressions indicates that the actual hfetime of a filament (as op- 
posed to the turbulent diffusion time across the edge layer) is virtually unbounded, provided 
direct damping due to resistivity or coUisional diffusion is negligible. During a filament life- 
time turbulence must diffuse across the scale a, many Ar-layer widths from the filament edge 
of to its center. However, in just a layer time r„ or rya^, the turbulence is reduced by many 
factors of e~^, while the filament density or current inside of the layer remains untouched. 
Consequently the width of the mixing layer at the edge of the filament continuously decreases 
even as the time to mix across it increases. The result is that mixing never extends to the 
filament core. This analysis shows that structures identified in the simulations as current 
filaments correlate spatially with a coherent density field, provided the density component 
is not destroyed by strong coUisional diffusion. 



6. Geometric and Statistical Properties 

The above analysis treats the current of the filament as locahzed. The current is max- 
imum at r = and becomes zero at r = a. This makes the shear of the filament magnetic 
field largest in the filament edge, and zero in the center. If it is true that the shear of this 
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field refracts turbulent KAW activity as described above, turbulence is suppressed where the 
shear is large. These properties are incorporate d in the spatial variation of a single quantity 
known as the Gaussian curvature (ITerryl 120001 ) . The Gaussian curvature is a property of 
vector fields that quantifies the difference between shear stresses and rotational behavior. In 
rectil inear coordinates the Gaussian curvature Ct of a vector field A(x, ?/) is (jMcWilliams 

igsj 

L'T = 7\ — + 1 7\ — ~ 7\ — • (,4oj 



For the total magnetic field in our cylindrical system this can be written 



dr\r 



r 36 



. dr\ r J r 09 



Jo+J 



(44) 



Inside the filament, turbulence is suppressed, and Ct is dominated by the filament field 
components Bg and Jq. Near the center, Jq is maximum and d{Bg/r)/dr = j' vanishes, 
making Ct negative. Toward the filament edge, j' becomes maximum as Jq goes to zero, 
making Ct positive. Outside the filament Ct is governed by bg, 6,., and j. These components 
must be roughly in balance. If they are not, the conditions for forming a coherent filament 
are repeated, and a structure should be present. Therefore, in regions where there are 
coherent filaments, the Gaussian curvature should have a strongly negative core surrounded 
by a strongly positive edge. Where there are no coherent structures the Gaussian curvature 
should be small. If this property is observed in simulations, it confirms the hypothesis that 
shear in the filament field refracts turbulent KAW activity in such a way as to suppress 
turbulent mixing of the structure. We note that the negative-core/positive-edge structure 
is predicted for current filaments of either sign, positive or n egative. This type of Gaussian 
curvature structure has been observed in recent simulations fjSmith fc Terryll2006l ). 



If the current filaments are well separated, their slow evolution relative to the decaying 
turbulence that surrounds them leads to a highly non Gaussian PDF. Assuming an initial 
PDF that is Gaussian with variance (J^), 



P(J) 



1 



/2^(J^) 



1/2 



exp 



-J2 



(45) 



it is possible to model subsequent evolution on the basis of the time scales derived pre- 
viously and the condition for strong refraction, Eq. (37). This condition stipulates that 
structures form where refraction is large, i.e., where mj' ^ i^lrya/d^- Since / ~ Jo/a and 
'0|r>a/a^ ~ {JaV^'^, structures occur for Jq > Jc = C(J^)^/^, where C is the smallest nu- 
merical factor above unity to guarantee strong refraction and mixing suppression. Given the 
latter, filaments reside on the tail of the PDF with high J and low probability P{J)dJ. 
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This probability is equal to the filament packing fraction, i.e., the fraction of 2D space occu- 
pied by current filaments. If, for simplicity, we assume that all filaments are of radius a, the 
faction of 2D space they occupy is {a/iy, where I is the mean distance between filaments. 
Therefore, 



r-Jc roo POO 2 

/ P{J)dJ+ P{J)dJ^2 P{J)dJ^(-), 

J -00 Jjc J Jc ^ ^ ^ 



(46) 



where we assume that P{J) is an even function. This expression gives the packing fraction 
as a function of the critical current Jc for filament formation. 

It is now straightforward to construct a heuristic model for the evolution from the initial 
distribution. The model applies for times that arc larger than the turbulent Alfvcn time, 
but shorter than the mean time between filament mergers. (Once filaments begin merging, 
their number and probability begin decreasing.) Prior to that time the filament part of the 
distribution with J > Jc is essentially unchanged, apart from the minor effects of slow erosion 
at the edge of the filaments. The probability that a fiuctuation is not a filament also remains 
fixed, but these fiuctuations decay in time. This means that the variance decreases while 
the probability remains fixed. The rate of decay is the turbulent Alfven time ta- Therefore 
the distribution can be written 



P(J,t) 



N{t) 



27r(J2)V2 
P(J,i) = 



exp 



-J2 



2(J2)expH/TA], 



2\l/2 



exp 



-J2 



for J < Jr 



for J > Jr 



(47) 



where ( J^) remains the initial variance, and N{t) is a time-dependent normalization constant 
that maintains J^" P(J, t)dJ at its initial value, i.e.. 



N{t) 





exp 






" dJ exp 


"- J2exp[VTA]/2(J2)" 



(48) 



The distribution P{J,t) becomes highly non Gaussian as t ^ because one part of the 
distribution (for J < Jc ) collapses onto the J = axis and becomes a delta function 5 (J), 
while the other part remains fixed. 



A simple measure of the deviation from a Gaussian distribution is the kurtosis, 

3jJ^P{J,t)dJ 



Jj^P{J,t)dJ 



2- 



(49) 
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The evolving kurtosis can be calculated directly from Eq. (47) . While the exact expression is 
not difficult to obtain, its asymptote is more revealing. The kurtosis diverges from the initial 
Gaussian value of 3 as the contribution from turbulent kinetic Alfven waves {J < Jc) decays 
and collapses to 6 (J). After a few Aflvcn times the kurtosis is dominated by the part with 
J > Jc, which, because it is stationary, represents the time asymptotic value for ta < t < tm- 
The time tm is the mean time to the first filament mergers. The time asymptotic kurtosis is 



to yield {£) exp{-J^/2{J^)) = (a//)Ve[l + 0((j2)3/2/j3)]_ xhe time-asymptotic kurtosis 



is much greater than the initial Gaussian value of 3 and is characterized by the initial value 
of the inverse packing fraction. Once filament mergers begin, the inverse packing fraction 
increases above the initial value (l/a)^. If [l(t)/a]'^ is the inverse packing fraction for t > tm, 
the above analysis suggests that the kurtosis will continue increasing as (3/2)[Z(i)/a]^ for 



We now consider the distribution of density. As shown in the previous section, the 
density present in the current filament also has suppressed mixing and is therefore coherent, 
or long lived. However, it is not spatially intermittent to the same degree as the current. 
Alfvenic dynamics indicate that n ~ -B, while Ampere's law stipulates that the magnetic field 
of the filament extends into the region r > a, falling off as r~^. Hence the density associated 
with filaments also is expected to fall off as for r > a. This spatially extended structure 
makes density less isolated than current. It produces higher probabilities for low values of 
density than those of decaying turbulence. This will yield a kurtosis closer to the Gaussian 
value of 3 than the kurtosis of the current. However, the distribution of low-level density 
associated with the structure likely will not be Gaussian. It is ultimately the distribution 
that matters for the scattering of rf pulsar signals. 

To construct the density PDF we seek the mapping of density onto the spatial area it 
occupies. We obtain this mapping for the filament density, assuming that the density of 
turbulence is low, and has effectively collapsed onto n = after a few ta, just as the current. 

The density for r > a goes as n = ano/r, where Uq is the value of the density at r = a. As 
shown in Fig. 3, the area occupied for a given density is 2'Krdr. This area is the probability 
when properly normalized, hence. 



Writing rdr in terms of ndn using n — ano/r, P{n) — Cnjr? , where C„ is the normalization 
constant chosen so that the probability integrated over the whole filament with its r~^- 




(50) 



late times. 



P{n)dn — 2Trrdr 



(51) 
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mantle equals the packing fraction, or probability of finding the filament in some sample 
area. With the long, slowly decaying tail of it is necessary to impose a cutoff to keep 
the PDF integrable. The cutoff, which we will label Uc, corresponds to the low level of 
decaying turbulence, but otherwise need not be specified. Consequently, the normalization 
is determined by 



for Tc < I, 



2Cr 



I 

J n. 



for Tc > I, 



(52) 



where rc, the radius at which n = ric, is Tc = ano/ric- The first of the two possibilities in 
Eq. (52) allows for a cutoff radius that is smaller than the mean distance between structures 
(of radius Vc), yielding a probability that is less than unity. If the cutoff radius is equal to 
or larger than the mean separation, then the structures are space filling and the probability 
is unity. Solving for C„, the normalized density PDF is 



P{n) 



or 



a Tie 
I no 



P{n) 



2^2 



rc>l 



or 



a Tic 
-> — 
I no 



(53) 



This distribution is defined for Uc < n < Uq. It captures only the contribution of filaments, 
and ignores the density inside r = a, which makes a small contribution to the PDF. 

This distribution is certainly non Gaussian, because it has a tail that decays slowly as 
n~3. However, depending on the length of the tail, which is set by n,. and no, the distribution 
may or may not deviate from a Gaussian is a significant way over Uc < n < no- This is 
quantified by the kurtosis. 



K{no,nc) = 



6Cn£'n\n-'dn) 



CI 



2/;° n2(n-3dn) 



2 ■ 



(54) 



Substituting from Eq. (53), we find that 



K{no, nc) 



3f 


1 - 


2nl/nl + ni/n^ 


4a2 


In (no/nc)J 



or 



a nc 
I no 



(55) 
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3nl 



1 - 2^2/^2 + ni/n^Q 



In {no/ Tic 



rc>l 



or 



a ric 
- > —. 
I no 



(56) 



These expressions are smaller than the current kurtosis by a factor 2[ln{no / nc)]'^ . Unless 
no/nc is quite large, the kurtosis may not rise significantly above 3. This is particularly 
true in simulations with limited resolution where dissipation will affect the density, either 
directly through a coUisional diffusion, or indirectly by resistive diffusion of current filaments. 
Kurtosis increases if nc decreases. However, while nc is tied to the decreasing turbulence level, 
regeneration of the turbulence by the mantle may prevent Uc from becoming very small. 
Nonetheless, mergers of filaments will decrease the packing fraction. Even if the density is 
space filling initially and satisfies Eq. (56), the mean filament separation will increase above 
Tc at some point, and the kurtosis will be given by Eq. (55). Then as the inverse packing 
fraction increases above P/a^, the kurtosis will rise. 

The n~^ falloff of the density PDF has intriguing implications for r/ scattering of pulsar 
signals. Noting that the scattering is produced by gradients of density, the extended density 
structure for r > a yields Vn = n' ~ l/r^. We can construct the PDF for n' following 
the procedure used for the PDF of n. Writing rdr in terms of n'dn' using r ~ (n')"^^"^, we 



recover 



P(n') 



\n 



/^2 ' 



(57) 



where Cn' is a constant. This is a Levy distribution , the type of distribution inferred in the 
scaling of pulsar signals (IBoldyrev fc Gwinnll2003bl ). Further exploration of the implications 
of these results to r/ scattering of pulsar signals remains an important question for future 
work. 



7. Conclusions 

We have examined the formation of coherent structures in decaying kinetic Alfven wave 
turbulence to determine if there is a dynamical mechanism in interstellar turbulence that 
leads to a non Gaussian PDF in the electron density. Such a PDF has been inferred from 
scalings in pulsar scintillation measurements. We use a model for kinetic Alfven wave tur- 
bulence that is applicable when there is a strong mean magnetic field. The nonlinearities 
couple density and magnetic field in the plane perpendicular the mean field in a way that 
is analogous to the coupling of fiow and magnetic field in reduced MHD. The model applies 
at scales on the order of the ion gyroradius and smaller. We show that the coherent current 
filaments previously observed to emerge from a Gaussian distribution in simulations of this 
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model (ICraddock et al.lll99ll ) result from strong refraction of turbulent kinetic Alfven waves. 
The refraction occurs in the edge of intense, localized current fluctuations, and is caused by 
the strongly sheared magnetic field associated with the current. This refraction localizes 
turbulent wave activity to the extreme edge of the filament, and impedes mixing (turbulent 
diffusion) of the filament current by the turbulence. From this analysis we conclud e that 
the turbulence suppression by sheared flows common in fusion plasmas ( Terryl 2000 ) has a 



magnetic analog in situations where there is no flow. This leads to a further conclusion that 
intermittent turbulence, which is generally associated with flows, can occur in situations 
where there is no flow. (By flow we mean ion motion. Electron motion is incorporated in 
the current.) We have derived a condition for the strength of magnetic shear required to 
produce the strong refraction and suppress mixing. We show that this condition yields a 
prediction for the Gaussian curvature of the magnetic field. This quantity is predicted to 
have large values inside the coherent current filaments, and small values everywhere else. 
Inside filaments the Gaussian curvature is negative at the center, and positive at the edge. 

The analysis shows that long-lived fluctuation structures form in the density and mag- 
netic field, provided damping is negligible. Like the current filaments, these structures avoid 
mixing because of the refraction of turbulent kinetic Alfven wave activity. Hence they oc- 
cur in the same physical location as the current filaments. However the localized nature 
of the current filaments give the long-lived magnetic field an extended region external to 
the current. In this region the field falls off as r^^, where r is distance from the center 
of the filament. Because kinetic Alfven wave dynamics yields an equipartition of density 
and magnetic field fluctuations, we posit that the long-lived density has a similar extended 
structure. As a result, the connection between coherent structure and localization that is 
true for the current, and makes it highly intermittent, does not apply to the density. While 
there is coherent long-lived density, it need not be localized. A sim ilar situation holds for 



vorticity and flow in 2D Navier Stokes turbulence (jMcWilliamslll984j ) . To explore this matter 
we have used the physics of the coherent structure formation to derive heuristic probability 
distribution functions for the current and density. As the turbulence decays, leaving intense 
current fluctuations as coherent current filaments, the kurtosis of current increases to a value 
proportional to the packing fraction. The kurtosis of density does not become as large, and 
could, under appropriate circumstances, remain close to the Gaussian value of 3. However 
mergers of structures in a situation with very weak dissipation could increase the kurtosis 
well above 3. More importantly, however, the density PDF is non Gaussian even when its 
kurtosis is not greatly different from 3. The structure external to the current gives the 
PDF a tail that goes as n~^. Mapping the structure to a PDF in density gradient, the 
density-gradient PDF decays as 1/n'^, a Levy distribution. This suggests that the mechanism 
described here may play a role in the scaling of pulsar r/ signals. 
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Several aspects of this problem need additional study. It is important to adapt these 
results to a steady state. Generally speaking, there is a dynamical link between decaying 
turbulence and turbulence in a stationary dissipation range. Hence these results, at least 
qualitatively, are relevant to the dissipation range. Dissipation begins to affect the spectrum 
at a scale th at is somewhat larger (order of magnitude) than the nominal dissipation scale 
(lFrischlll995l ). Structures such as these would correspond to active, filamentary regions of 
dissipation analogous to those observed in neutral gas clouds, albeit at a much smaller scale 
and with no accompanying flow shear signature. Intermittent structures can extend into 
the stationary inertial range, but the analysis presented here must be modified. In the 
inertial range, turbulence is replenished, allowing the slow mixing of a coherent structure 
to continue until it is gone. Structures are also regenerated by the turbulence, and the 
statistics is ultimately set by a balance of mixing and regeneration rates. The mixing rates 
calculated here are sufficiently slow in strong filaments, that coherent structure formation 
is expected even in a steady state. There is also a possible link between structures in the 
larger scale range of shear Alven excitations and KAW excitations. These questions will be 
explored in future work. While gyr or adius- scale KAW turbulence may arise in astrophysical 
contexts other than the ISM, the small scales make it unlikely that astrophysical observations 
will be available for testing this theory. Therefore, simulations should be used to check key 
conclusions from the theoretical work presented here. These i nclude the formation of density 
structures, which was not reported in ICraddock et al.l (Il99ll ). the structure of the Gaussian 
curvature, which validates the refraction hypothesis, and the existence of the structure in 
the density and its effect on the PDF. The effect of this type of density field on r/ scattering 
remains the underlying question, and modeling of the scattering with simulated fields should 
be pursued. 
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A. Appendix 



Closures truncate the moment hierarchy that is generated when averages are taken of 
nonlinear equations. The closure we have used is of the eddy damped quasi normal Mark ovian 
variety, and follows the steps of the closure calculation described in (ITerry et al.ll200ll ). The 
nonlinear decorrelation is calculated consistent with the statistical ansatz, not imposed ad 
hoc. The closure equations are given in Eqs. (23) and (24). The other diffusivities not given 



in 



Eq. (25) 



are 
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^3("^, 7) = 7)^2(^,7) + - 7 + L'^v("^'7)^V,^„ - ——noj , (A7) 

and Pm and Ki{m,'j) are given in Eqs. (26) and (27). These expressions contain both hnear 
wave terms and nonhnear diffusivities, and are vahd in both weak and strong turbulence 
regimes. Outside filaments, where turbulence levels are evaluated to derive the strong re- 
fraction condition, Eq. (37), the turbulence is strong. The strong turbulence limit of the 
above expressions yields the diffusivity d that appears in Eq. (37). 
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Fig. 1.— 



Distortion of a kinetic Alfven wave phase front by a slieared filament field. 
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Fig. 2. — Boundary layer at edge of filament (just inside r — a). External fluctuations 
rapidly decay across the layer. 
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Fig. 3. — The probability of density outside r = a maps onto the annular area 27rrdr. 



